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We consider energetics and structural properties of a many particle system in one dimension with 
pairwise contact interactions confined in a parabolic external potential. To render the problem 
analytically solvable, we use the harmonic approximation scheme at the level of the Hamiltonian. 
We investigate the scaling with particle number of the ground state energies for systems consisting 
of identical bosons or fermions. We then proceed to focus on bosonic systems and make a detailed 
comparison to known exact results in the absence of the parabolic external trap for three-body 
systems. We also consider the thermodynamics of the harmonic model which turns out to be similar 
for bosons and fermions due to the lack of degeneracy in one dimension. 

PACS numbers: 



I. INTRODUCTION 

Exactly solvable models are key players in few- and 
many-body quantum mechanics and are indeed also fas- 
cinating creatures 0, 0|. The analytical intractability 
of general iV-body systems makes it extremely impor- 
tant to have exact results to allow benchmark tests of 
complicated numerical methods. Unfortunately, exactly 
solvable models are few and far between and are more of- 
ten found in low-dimensional systems. In the case of one 
spatial dimension, the famous Bethe ansatz Q was suc- 
ccsfully applied to models of bosons with zero-range in- 
teractions [4 6] , and later to interactions with long range 
0,1- 

Meanwhile, an exciting direction in the field of cold 
atomic gases aims at the study of low-dimensional sys- 
tems in general and one-dimensional setups in particu- 
lar [9j. A noticable highlight of this pursuit is the ex- 
perimental realization of the so-called Tonks-Girardeau 
gas [l0|-[llj], where strongly-interacting one-dimensional 
bosons become impenetrable objects and behave similar 
to fermions [HI, [lij . Very recently, it has even become 
possible to study this interesting regime in the limit of 
small particle number [IH, [l6[ . 

Here we study an exactly solvable model of an TV-body 
system in an external parabolic confinement. The outer 
trap is always present in cold atomic gas experiments 
but can often be neglected or treated in a local density 
approximation when large systems are studied. However, 
for smaller particle numbers the effect of the outer trap 
becomes important for the structure and dynamics of the 
system. To make the iV-body problem tractable we use a 
harmonic Hamiltonian approximation [§} with carefully 
chosen parameters that reproduce essential features of 
two atoms interacting via short-range interactions in a 
parabolic trap [l7j . 

The paper is organized as follows. After a discussion of 
the harmonic methods we provide details on how the pa- 
rameters of the harmonic interactions are obtained from 
knowledge of the exact solution of the two-body problem 
originally obtained by Busch et al. [l8| . We then present 
results for the ground-state energies and the radii. Our 



focus is on bosonic systems, but we do present a few re- 
sults for fermions as well. A comparison to the exact re- 
sults of MacGuire @ is then made in relevant limits. We 
also compute the one-body density matrix and its largest 
eigenvalue to obtain the condensate fraction at zero tem- 
perature. Lastly, we discuss the thermodynamics of our 
model and then proceed to conclusion and outlook for 
future work. 



II. METHOD 

We consider a system of N quantum particles interact- 
ing pairwise via a delta function interaction and confined 
by a harmonic external potential in one-spatial dimen- 
sion. The Hamiltonian for this system is 
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where m is the mass of the particles, ujq is the frequency 
of the external field, and a is the one-dimensional scatter- 
ing length which parameterizes the strength of the two- 
body interaction (we will discuss its relation to the three- 
dimensional scattering length below). The external field, 
wo i defines the length scale of our system, I 2 = hi (mwp). 
For two particles, ([I]) has been solved in, e.g., jla. fl9|. 
and their results of two-body energies and wave functions 
for a given scattering length are used to determine the 
parameters of our model. We consider only the bound 
molecular branch of the system, i.e., a > 0. 

In our general harmonic approximation scheme we re- 
place Hamiltonian ([TJ with 
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where /i = m/2 is the reduced mass of the two-body 
system, LOi n is the interacting frequency, and Vs is an 
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energy shift. The solution to this equation for N particles 
for general systems is described in detail in 1171 and more 
specifically for identical particles in H and [2dj|. 

The parameters Ui n and Vs are now chosen to fit per- 
tinent properties of the original Hamiltonian at the two- 
body level in the usual spirit of constructive descriptions 
of iV-body systems based on two-body interactions. Here 
we impose the condition that at the two-body level the 
harmonic oscillator reproduces the energy and average 
square radius of the exact two-body solution. In order to 
fulfil these constraints, we use the size of the system to 
determine the interaction frequency through the relation 



(3) 



where x = X\ — X2 is the relative coordinate in the two- 
body system. The energy shift is determined by requiring 
that the model reproduces the energy of the two-body 
system 



E2 = 2 H \J"in+"o + Vs- 
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The quantities E2 and (x 2 ) can be easily obtained by nu- 
merically solving the transcendental equations fulfilled 
by the exact solution for two bosons (or two fermions in 
different spin states) interacting via a zero-range interac- 
tion in a parabolic trap [l8| . Note that when we consider 
identical fermions below, there can in principle only be 
a non-zero two-body interaction in states that are odd 
under exchange of the two fermions. In this paper we are 
mostly interested in the scaling behaviour with particle 
number of the fermions as the interaction frequency is 
varied (through the scattering length) and we therefore 
ignore this point and use the same input for fermions as 
for bosons. The most important difference is of course 
the quantum statistics which is taken fully into account. 

We note that in cold atoms, effective one-dimensional 
setups are created by using a tightly confining potential 
in two transverse directions, usually through the applica- 
tion of an optical lattice 0. In a deep transverse lattice, 
the atoms are then effectively only occupying the lowest 
transverse motional degree of freedom. This transverse 
degree of freedom, however, has an influence on the effec- 
tive one-dimensional scattering length that described the 
inter-atomic interaction in the system. A mapping exists 
between the true three-dimensional scattering length and 
the effective one-dimensional scattering length that takes 
the transverse degrees of freedom explicitly into account 
[2lj . This mapping should thus be applied before com- 
parison to realistic experiments. 

Once the parameters have been determined, ((2} can be 
solved. The ground state energy for identical bosons is 

Ef) = \N{N - 1)V S + ±(N - l)hu r + ±fk* , (5) 



where 



is the (N— l)-degenerate frequency that comes out of the 
solution of © [17j • The last term in (O is the energy re- 
lated to the center of mass motion and will be ignored as 
we are interested in the internal dynamics of the system. 

Though most of this work concerns bosons, input for 
identical fermions is also possible. Since fermions obey 
the Pauli Principle, N — 1 particles must be placed in 
higher oscillator levels as there are no degeneracies in 
one dimension. We can obtain the amount of energy in 
the ground state by filling consecutive oscillator levels, 
i.e. 
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The complete ground state energy is then 
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This is a different scaling than for the bosons in ([5])- 
In the fermion case, the energy is dominated by the 
N 2 huj r /2 term which scales as N 5 / 2 . 
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FIG. 1: Top panel: Energies per particle for bosons confined 
on one dimension. The j/-axis is stopped at -20 in order to 
show the general behaviour of many curves, though the short- 
est scattering lengths quickly leave the figure. The scattering 
lengths plotted, moving from bottom to top are 0.4., 0.5, 1.0, 
2.0, 5.0, 10.0, and 100.0 Bottom panel: Energies divided by 
N 2 for the bosons to show their scaling behaviour with N. 
The order of the scattering lengths from bottom to is the 
same as in the top panel. 



III. RESULTS 

We now present numerical results for systems with 
N = 3 — 30 particles, studying their energetics, the ra- 
dial behaviour of the systems, and the one-body density 
matrix. Here we focus on the case of identical bosons. 
A particular issue is the ratio of the three-body energy 
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to the two-body energy. An exact formula of MacGuire 
@ applies to one-dimensional systems of bosons with 
zero-range interactions and we make a comparison of the 
harmonic results to that model in the strongly-bound 
limit where the two-body energy is large and negative. 
Lastly, we consider also the thermodynamics of the sys- 
tem within the harmonic approximation. 



A. Energies and radii 

The ground state energy for bosons is shown in figure 
Q] and that for fermions in figure [2] for different particle 
numbers and scattering legnths. Notice that the upper 
panels show the energy per particle, while in the lower 
panels the energy is divided by a different power that 
we will discuss below. A striking feature to notice is the 
positivity of the energy for fermions, while that of the 
bosons can have both signs. This is a consequence of the 
Pauli principle which implies that the fermions will have 
to occupy higher orbitals for lack of degeneracies in one 
dimension as discussed above. Even in the case of very 
small scattering lengths (and thus large and negative two- 
body binding energies) the contribution from the higher 
orbits makes the overall ground state energy positive. 

Once the two-body problem is solved and the param- 
eters are determined, the ground state energy for bosons 
is given by (JS|) . As seen in figure [TJ the sign of the en- 
ergy of the Bose system is mostly negative, but for small 
systems at large scattering lengths, the energy turns pos- 
itive. This is understood from the fact that the energy 
shift, Vs, is negative and the shift term in the energy of 
(O scales with A 2 while the positive oscillator frequency 
term scales with the lower power A 3 / 2 . For bosons, the 
energy is positive for all particle numbers examined for 
a 1 I = 100. Using ([3]), one can derive the critical number 
for bosons where the energy changes sign. It is given by 

ft 2 u4 + 2V^a4/4 + 4V^) 

Ncrit - W i . (9) 

From this relation one can calculate that for bosons the 
energy becomes negative at a/l = 100 when A > 116. 

The bottom panels of Figures Q] and [2] show the energies 
divided by the relevant scaling factor (A 2 for bosons and 
A 5 / 2 for fermions) to emphasize the scaling behaviour 
at large A. In the figures we can see that the asymp- 
totic scalings are nicely approached when the scattering 
length is not too small. For the case of small scatter- 
ing length, the two-body energy will in general behave as 
E 2 ~ — 1/a 2 while (x 2 ) ~ a 2 , since in this limit the trap 
can be ignored on the bound state branch in the spec- 
trum that we study here. However, the relations in ([3]) 
and ((H) that we use to determine our oscillator param- 
eters now imply that both Vs and LUi n must scale with 
a~ 2 to be fulfilled (remember that E 2 is negative). This 
means that there will be a competition between the terms 
in the ground state energies given by (J5J) and ([5]). The 
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FIG. 2: Top panel: Energies per particle for fermions confined 
on one dimension. The plotted scattering lengths are in the 
opposite order from the bosons as moves from the bottom to 
the top of the plot. Bottom panel: Fermion energies divided 
iV 5 / 2 to clearly show the scaling of the energy with N. 



asymptotic behaviour for large A is therefore approached 
more slowly for small a. 

We also calculate the relative size of the bosonic ground 
state wave functions, ({X — Xcm) 2 ), as a function of 
particle number for several different scattering lengths. 
Note here that X denotes the single particle coordinate 
of one of the bosons and that no index is needed since 
the particles are identical. The results are shown in figure 
[U One can see for the large scattering lengths that the 
radius increases before eventually decreasing for larger 
particle numbers (more than 20 for a/l — 100). For 
scattering lengths of two and smaller, the size decreases 
monotonically with increasing particle number. Again 
this is connected to the fact that the small a regime has 
very strong two-body binding and thus small (x 2 ) which 
is imprinted on the A-body system which tends to be 
very compact and presumably leads to strong clusteriza- 
tion in the real system followed by loss of atoms from the 
trap. 

For fermions (not shown), the Pauli principle and the 
subsequent need to occupy higher orbitals means that the 
system will in general be larger for the same scattering 
length and will not show a decreasing behaviour as A is 
increased, but rather a slight increase. This can be seen 
from the energetics of the ground state in ([5]) which is 
dominated by the positive term containing the oscillator 
frequency for large A. This implies that the radius will 
stay large, constrained solely by the external trap. 

We notice that the dimensional effect on the radius is 
quite clear when comparing to previous studies of bosonic 
systems in two and three dimensions [l7j . In the current 
one dimensional study, we see a non-monotonic behaviour 
in the radii as function of A which can also be seen in 
two dimensions but which is almost completely absent in 
the corresponding three-dimensional system. The scat- 
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FIG. 3: Radius (x — Xcm) 2 for bosons as a function of 
particle number for several different scattering lengths. The 
scattering lengths are, from bottom to top, 0.4, 0.5, 1, 2, 5, 
10, and 100. 



tering lengths for which the maximum in the size occurs 
are roughly those where the three-body energy comes out 
positive. This implies an increased radial size, although 
still restricted by the external trapping potential. For 
large N, however, the bosons will always become nega- 
tive in energy as the shift term dominates in ([5]) and the 
radius goes down again. The fact that lower dimensions 
show a more pronounced non-monotonicity of the size can 
then be traced to the fact that the positive energy con- 
tribution from zero-point motion grows with dimension 
and washes out the behaviour. For fermions, as discussed 
above, the energy is positive from the start and the radius 
stays large. However, the degeneracies allowed in higher 
dimensions adds extra ingredients and shell-structure to 
both energy and to the radial size. 



B. Three boson energies 

We now consider the case of three bosons within our 
harmonic approximation scheme in order to make a 
comparison to the results obtained by MacGuire [6[ in 
the absence of external confinement. The exact bound 
state energy, En, of N bosons interacting via attrac- 
tive pairwise zero-range interactions in a homogeneous 
one-dimensional space can be expressed in term of the 
two-body energy, E 2 , as 

E N = ^N(N 2 -1)E 2 . (10) 

For the case of three particles this becomes E3 = AE2- 
Since we have an external trap, we do not expect to re- 
produce this result. However, in the limit a — > where 
E-2 — > —00 the external trap should become negligible 
and a comparison can be made. 



We plot the ratio E3/E2 as a function of scattering 
length in figure 0] and as a function of E 2 in figure [5] 
The plot ends at a/l ~ 0.25 at which point it becomes 
numerically very challenging to compute the wave func- 
tion based on the exact solution of Busch et al. (l8| . 
However, beyond that we can use the exact solution as 
discussed below. 

One can see in figure 2] as the scattering length de- 
creases (more clearly seen as E2 — > — 00 in figure [5]), that 
the ratio approach a limit. Using ([5]) and (U]) in the limit 
where a — > and io r — > uj in (see ©), we get that 

t-'+G/i-§)£- 

Relating uj in to (x 2 ) through ([3]), we obtain 

<i2 » 

We can now use the exact wave function for a delta func- 
tion potential in one dimension, ip{x) = Ae~ KX where 
k = ^-2^E 2 /h 2 , to obtain ti 2 /(2fiE 2 (x 2 )) = 2. Our 
limit thus becomes E3/E2 ~ 3.55. This is also the num- 
ber we find numerically for the smallest value of a that we 
could access. Comparing to the exact value of MacGuire, 
E3/E2 = 4, this implies that the energetics of our model 
is accurate to around 10% for small scattering lengths. 
The wave function in the harmonic model is a gaussian 
which in the limit a — >■ will tend to a delta function, 
similar to the exact solution, so we also expect the har- 
monic model to provide an accurate structural descrip- 
tion in this limit. 

On the other hand, when a becomes very large (the 
unitarity limit) the trap plays an important role. Looking 
back at the basic Hamiltonian in ([I]), we see that when 
a — > 00, the interaction term tends to zero and we should 
be dominated by the trap only. However, the exact result 
does show that there is a shift of the energy in this limit so 
that E2 = huj/2 for a — > 00 (the non- interacting system 
has twice the energy since there is zero-point motion from 
both particles). Within our model we have uj r — > ujq and 
(x 2 ) — > l 2 /2 when a — > 00 as the external field provides 
the only length scale left in the problem. We have 

^=3-^, (13) 
E2 2E2 

and in the large a limit, E2 — > Hujq/2, so we obtain a 
limiting value of £3 — > 2. At a/l = 100, we numerically 
obtain 1.98 for this ratio so we reproduce this limit. The 
wave function becomes essentially exact in this limit, and 
we thus expect the structure to be well reproduces also. 

In the limit of large N, we can also estimate the scal- 
ings and compare to the exact results when a — > 0. From 
([5|) we have that the bosons scale as En ~ N(N — l)i?2, 
which means that we have an underbinding by a factor 
of N + 1. Fermions in ([5]) are different since here we have 
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FIG. 4: The ratio E3/E2 as a function of scattering length 
for bosons. We plot only the cases where both E3 and E2 are 
negative. 



FIG. 6: The value of A at several values of the ID scattering 
length. From bottom to top they are 0.4, 0.5, 1, 2, 5, 10, and 
100. 



UJ 
LU 



6 
5.5 

5 
4.5 

4 
3.5 

3 



1 1 1 — 

16 -14 -12 -10 



C. One-body density 



-8 -6 

E 2 (hco ) 



FIG. 5: The ratio E3/E2 as a function of E2 for bosons. We 
plot only the cases where both E3 and E2 are negative. 



En ~ N 5 ^ 2 E 2 , i.e. the underbinding is only by a fac- 
tor of about %/iV. Notice, however, that in the opposite 
limit of a — > 00 , the dominant term in ([5]) and © change 
as Vs — > and ui r — > ujq. This is an essentially non- 
interacting situation and we have the intuitively obvious 
result E N ~ Nhuj /2, which is just the scaling dictated 
by the trap. For values of a away from these two lim- 
its, the scaling is an interplay of both terms and we get 
a more complicated behaviour. We note that one could 
also turn these arguments upside down and use the ex- 
act results for TV particles in the a — > limit to fit the 
parameters of the model, for instance by using E 3 = AE 2 
instead of either the frequence of radius condition we im- 
pose on the two-body problem. Overall, we expect this to 
provide only minor quantitative changes in the harmonic 
model and its predictions. 



The one body density matrix 



p{xi,x' 1 ) — N I tp*(xi,.. .,x N )ip(x 1 , . . .,x N )dx 2 ■ ..dx N , 



(14) 

is the most basic correlation function and is the starting 
point of statistical property calculations. It contains in- 
formation about the mean-field nature of the wave func- 
tion for a bosonic system. It is of interest in itself, but we 
will focus on its largest eigenvalue, A, which directly mea- 
sures how much a given state has the structure of a coher- 
ent state. Figure [6] shows A as a function of the number 
of bosons for several scattering lengths. We mainly see 
an increase with particle number, though there are slight 
decreases for small particle numbers and large scattering 
lengths. In contrast with our previous results in 17[ for 
higher dimensions, only scattering lengths of order unity 
or smaller show a A much different than one, showing 
that the persistence of the mean-field structure provided 
by the external potential is quite strong in one dimen- 
sion implying that the mean-field approximation is very 
good. The difference to higher dimensions arise from the 
degeneracies and from the reduction implied by the fact 
that one can factorize A, so that for instance A3/5 = A 3 
(see [13). 



D. Thermodynamics 

The energy spectrum can be obtained by considering 
center-of-mass excitations in the system and also the in- 
ternal excitations characterized by the frequency in ([B]). 
To address the thermodynamics of the systems, one can 
proceed as described in detail in [20|. One naturally 
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starts with the partition function 



Z(T) = J29iexp[-Ei/(k B T)], 



(15) 



where E^ is the energy of the ith state, gi is the de- 
generacy of that state, and ks is Boltzmann's constant. 
ID systems are simpler since there is no degeneracy in 
the energy levels in contrast to higher dimensions. With 
the sequence of energies and single-particle degenera- 
cies, thermodynamic quantities can be calculated. Once 
the spectrum and degeneracies of our systems have been 
worked out, we can calculate basic thermodynamic quan- 
tities in the canonical ensemble. In ID, the number of 
states grows much more slowly than in higher dimen- 
sions, though one does have to climb somewhat high in 
the spectrum due to the lack of degeneracy. In fact, in 
ID there is no difference between bosons and fermions in 
the sequence of the degeneracies of the A-body excited 
states. This comes from the fact that all the Pauli princi- 
ple implies for fermions is that the ground state has them 
occupying higher orbit which increases the ground state 
energy in comparison to bosons which can all go into the 
lowest orbital. Excitations on top of these ground states 
will now cost the same amount of energy and they will 
all have the same degeneracy since it is simply a matter 
of distribution of a number of particles in levels to obtain 
a specific total energy cost (the energy of the excitation 
above the ground state). 

This latter fact implies that some thermodynamic 
quantities such as the entropy and heat capacity are 
identical for bosons and fermions at the same scatter- 
ing length. Working in the canonical ensemble, in figure 
[7] we show the heat capacity C, which is the same for 
bosons and fermions. The heat capacity is related to the 
partition function by 
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(16) 
(17) 



For all particle numbers, we see an initial increase, which 
then slows down at higher temperatures. This initial in- 
crease is the saturation of the center of mass mode. The 
temperature, where the slope changes, depends on par- 
ticle number since more particles present in the system 
implies that the center of mass mode comprises a smaller 
share of the overall excitation energy. At the larger scat- 
tering lengths, the heat capacity then gradually increases 
towards the equipartion limit value of one. At small scat- 
tering lengths, the large size of the excitation frequency 
sharply delays the approach to the cquipartition limit. 
For the same scattering length, more particles mean a 
slower approach to the high temperature limit. This is 
clear from ©, which shows that the excitation frequency 
increases with particle number. 
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FIG. 7: The heat capacity per particle at several boson num- 
bers and scattering lengths, identified as TV, a, as a function 
of temperature. 



IV. SUMMARY AND OUTLOOK 

We have presented results for an A-body system of 
particles moving in one dimension and interacting via 
pairwise zero-range interactions under the influence of 
an external harmonic oscillator trap. The Hamiltonian 
is solved analytically within the harmonic approximation 
where the parameters of the model are adjusted to repro- 
duce properties of the exact solution of the corresponding 
two-body problem. This is along the lines of the stan- 
dard approach to many-body problems in many fields of 
physics. 

The ground state energies for both bosonic and 
fermionic systems were determined and we discussed 
the scaling of these quantities with the particle number. 
We also obtained the radial size of the bosonic systems 
and compared to previous findings in two- and three- 
dimensional setups. A careful study of the energy of three 
bosons in the harmonic model reveals that in the limit of 
large two-body binding energy where the external trap 
can be neglected, the harmonic model results are about 
10% below the exact result abtained by MacGuire [y], 
while in the opposite limit of negligible two-body binding 
energy (large scattering length) the model becomes essen- 
tially exact. For larger particle numbers, the harmonic 
approximation tends to underbind the system compared 
to the exact results when the two-body system is strongly 
bound. 

Additionally, we studied the one-body density matrix. 
The main goal was the computation of its largest eigen- 
value to characterize the degree of coherence in the sys- 
tem. The results turn out to be similar to those ob- 
tained in two and three dimensions, except that the one- 
dimensional coherence turns out to be more robust except 
for large two-body binding energies. Lastly, we studied 
the thermodynamics of the one-dimensional system. The 
specific heat was found to behave in a manner that is also 
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quite similar to higher dimensions. However, it is worth 
noticing that within the harmonic model, the spectrum 
of excitation for bosons and fermions turns out to be ex- 
actly the same. Due to the lack of degeneracy of the 
levels in one dimensions, the only difference is therefore 
found in the ground state energy. 

The one-dimensional setup studied here provides a 
very valuable benchmark for the harmonic models that 
have been used in higher dimensions in different con- 
texts d, [23l-l25| . Another good test ground for the har- 
monic models is systems that have long-range interac- 
tion. In particular, systems of dipolar molecules in one- 
dimensional tubes have been studied recently due to their 
interesting few- and many-body structure [26l435l |. The 
dipole-dipole interaction in two different one-dimensional 
tubes is typically characterized by having a pocket which 
becomes deep when the dipolar interaction is strong [32j . 
This implies that a harmonic approximation is a good 
starting point. This approach has been studied in the 
case of two-dimensional layers with dipolar particles (36l — 



|39| and compares very well to exact results [40l - |42} away 
from the weak-coupling limit. The application of the 
harmonic approximation to one-dimensional arrays con- 
taining dipolar particles is a topic for future work. 

Another interesting direction in which the harmonic 
model can be used is for the study of higher-order terms 
in the interaction such as those coming from effective- 
range corrections to the model of Busch et al. [H, |43[ . 
Recent experiments have shown a potential need for such 
corrections in both two- [44[ and three-dimensional cold 
atomic gas systems [45|. A possible way to account for 
such corrections would be to include quartic interaction 
terms an then carefully fit the additional parameter to 
few-body properties. Within the harmonic model, the 
quartic interactions can be computed quite straightfor- 
wardly through gaussian integrals. This could then be 
compared to mean-field results that include effective- 
range corrections in the limit of large particle numbers 



[1] Mattis D 1993 The Many-Body Problem (Singapore: 

World Scientific) 
[2] Sutherland B 2004 Beautiful Models (Singapore: World 

Scientific) 
[3] Bethe H 1931 Z. Physik 71 205 
[4] Lieb E H and Liniger W 1963 Phys. Rev. 130 1605 
[5] Lieb E H 1963 Phys. Rev. 130 1616 
[6] MacGuire J 1964 J. Math. Phys. 5 622 
[7] Sutherland B 1971 J. Math. Phys. 12 246 
[8] Calogero F 1971 J. Math. Phys. 12 419 
[9] Bloch I, Dalibard J and Zwerger W 2008 Rev. Mod. Phys. 

80 885 

[10] Parades B, Widera A, Murg V, Mandel O, Foiling S, 
Cirac I, Shlyapnikov G V, Hansen T W and Bloch I 2004 
Nature 429 277 

[11] Kinoshita T, Wenger T and Weiss D S 2004 Science 305 
1125 

[12] Haller E, Gustavsson M, Mark M J, Danzl J G, Hart R, 

Pupillo G and Nagerl H-C 2009 Science 325 1224 
[13] Tonks L 1936 Phys. Rev. 50 955 
[14] Girardeau M 1960 J. Math. Phys. 1 516 
[15] Serwane F, Ziirn G, Lompe T, Ottenstein T B, Wenz A 

N and Jochim S 2011 Science 332 6027 
[16] Ziirn G, Serwane F, Lompe T, Wenz A N, Ries M G, 

Bohn J E and Jochim S 2012 Phys. Rev. Lett. 108 075303 
[17] Armstrong J R, Zinner N T, Fedorov D V and Jensen A 

S 2011 J. Phys. B 44 055303 
[18] Busch T, Englert B-G, Rzazewski K and Wilkens M 1998 

Found. Phys. 28 549 
[19] Farrell A and van Zyl B J. Phys. A 2010 43 015302 
[20] Armstrong J R, Zinner N T, Fedorov D V and Jensen A 

S 2012 Phys. Rev. E 85 021117 
[21] Olshanii M 1998 Phys. Rev. Lett. 81 938 
[22] Yang C N 1962 Rev. Mod. Phys. 34 694 
[23] Zaluska-Kotur M A, Gajda M, Orlowski A and 

Mostowski J 2000 Phys. Rev. A 61 033613 
[24] Yan J 2003 J. Stat. Phys. 113 623 
[25] Gajda M 2006 Phys. Rev. A 73 023603 



[26] Argiielles A and Santos L 2007 Phys. Rev. A 75 053613 
[27] Kollath C, Meyer J S and Giamarchi T 2008 Phys. Rev. 

Lett. 100 130403 
[28] Chang C-M, Shen W-C, Lai C-Y, Chen P and Wang D-W 

2009 Phys. Rev. A 79 053630 
[29] Huang Y-P and Wang D-W 2009 Phys. Rev. A 80 053610 
[30] Wunsch B, Zinner N T, Mekhov I B, Huang S-J, Wang 

D-W and Demler E 2011 Phys. Rev. Lett. 107 073201 
[31] Dalmonte M, Zoller P and Pupillo G 2011 Phys. Rev. 

Lett. 107 163202 
[32] Zinner, N T, Wunsch B, Mekhov I B, Huang S-J, Wang 

D-W and Demler E 2011 Phys. Rev. A 84 063606 
[33] Fellow J M and Carr S T 2011 Phys. Rev. A 84 051602(R) 
[34] Knap M, Berg E, Ganahl M and Demler E 2012 Phys. 

Rev. B 86 064501 
[35] Lecheminant P and Nonne H 2012 Phys. Rev. B 85 

195121 

[36] Wang D-W, Lukin M D and Demler E 2006 Phys. Rev. 

Lett. 97 180413 
[37] Armstrong J R, Zinner N T, Fedorov D V and Jensen A 

S 2010 Europhys. Lett. 91 16001 
[38] Zinner N T, Armstrong J R, Volosniev A G, Fedorov D 

V and Jensen A S 2012 Few-Body Syst. 53 369 
[39] Armstrong J R, Zinner, N T, Fedorov D V and Jensen A 

S 2012 Eur. Phys. J. D 66 85 
[40] Volosniev A G, Fedorov D V, Jensen A S and Zinner N 

T 2011 Phys. Rev. Lett. 106 250401 
[41] Volosniev A G, Zinner N T, Fedorov D V, Jensen A S 

and Wunsch B 2011 J. Phys. B 44 125301 
[42] Volosniev A G, Fedorov D V, Jensen A S and Zinner N 

T 2012 Phys. Rev. A 85 023609 
[43] Zinner N T 2012 J. Phys. A 45 205302 
[44] Baur S K, Frohlich B, Feld M, Vogt E, Perlot D, Koschor- 

reck M and Kohl M 2012 Phys. Rev. A 85 061604(R) 
[45] Hazlett E L, Zhang Y, Stites R W and O'Hara K M 2012 

Phys. Rev. Lett. 108 045304 
[46] Fu H, Wang Y and Gao B 2003 Phys. Rev. A 67 053612 
[47] Zinner N T and Th0gersen M 2009 Phys. Rev. A 80 



8 



023607 [49] Th0gersen M, Zinner N T and Jensen A S 2009 Phys. 

[48] Zinner N T 2009 Preprint larXiv:0909.'T3"T4l Rev. A 80 043625 



